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CN ; ABSTRACT 

■ 

. Aims. HESS J1745-303 is an extended, unidentified VHE (very high energy) gamma-ray source discovered using H.E.S.S. in the Galactic Plane 
qq ' Survey. Since no obvious counterpart has previously been found in longer-wavelength data, the processes that power the VHE emission are not 
well understood. 

Methods. Combining the latest VHE data with recent XMM-Newton observations and a variety of source catalogs and lower-energy survey data, 
we attempt to match (from an energetic and positional standpoint) the various parts of the emission of HESS J1745-303 with possible candidates. 
Results. Though no single counterpart is found to fully explain the VHE emission, we postulate that at least a fraction of the VHE source may be 
explained by a supernova-remnant/molecular-cloud association and/or a high-spin-down-flux pulsar. 
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1. Introduction 

The Galactic center region (roughly between |/| < 2°, \b\ < 1°) 
is densely populated with possible VHE emission candidates: 
supernova remnants (SNRs), dense molecular clouds, pulsar 
wind nebulae (PWNe), X-ray binaries (XRBs), and a variety 
of unidentified sources seen in lower wavebands. In this re- 
gion, VHE emission has been detected from an as yet uniden- 
tified point-like source at the Galactic center (possibly asso- 
ciated with the super- m assive black hole Sg r A* or a PWN) 
dTsuchiva et all 12004 iKosack et al.l 120041 lAha 
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(lAharonian et alj2005al) (point-like for H.E.S.S.), and from a re- 
gion of diffuse emission approximately ±1° in longitude, which 
is most likely associated wit h the interaction of cosm ic-ray parti- 
cles with molecular clouds (lAharonian et al ] l2006d) . and finally 
HESS J1745-303, an extended, unidentified VHE gamma-ray 
source lying approximately a half degree below the Galactic 
plane at / = -0.4. HESS J1745-30 3 was first discovered i n 
the H.E.S.S. Galactic Plane Survey ( lAharonian et all l2006dh . 
Subsequent observations of the region using H.E.S.S. have pro- 
vided increased exposure of this object, and thus a more detailed 
study is now possible. 

VHE gamma rays are typically thought to be produced via 
two general mechanisms: the up-scatter of lower-energy pho- 
tons by high-energy electrons via the inverse-Compton process, 
or the production and subsequent decay of 7r°s produced in the 
interactions of high-energy hadrons. Though it is difficult to dis- 
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tinguish between a purely leptonic or hadronic scenario in many 
of the currently published VHE sources, in cases where it is 
known that cosmic rays are interacting with a dense medium, 
the hadronic scenario becomes more viable. This is true for 
example on the Galactic center ridge, where the VHE emis- 
sion is see n to roughly follow the location of dense molecu- 
lar clouds (lAharonian et alJl2006d) . or possibly in the case of 
SNRs embe dded in dense regions o f the interstellar medium 
(ISM) (e.g. lAharonian et alj 2006d) . The flux of gamma rays 
produced via the hadronic production of 7r°s depends linearly 
on the density of the surrounding medium. For typical Galactic 
SNRs (few kpc distance, age of > 1000 yrs, in a medium with 
average density ~ 1 cm -3 ), the flux above 1 GeV is predicted to 
be quite small (below t he EGRET or even H.E.S.S. sensitivity) 
( lAharonian et al. I fl99l . However, in the case where the super- 
nova shock is interacting with a dense molecular cloud (as in the 
case of G359. 1 -0.5, discussed later), the emission can be signif- 
icantly enhanced. Moreover, the expected flux in the TeV energy 
range from an interacting SNR shock may be much higher than 
the extrapolation of the spe ctrum measured in the GeV range 
( lAharonian & Atovanll 19961) . 

Though the region ar ound HESS J1745- 303 is well cov- 
ered by radio (e.g. VLA. lLaRosa et al.ll2000l) . and X-ray (e.g. 
ROSAT. rVoges et aT]|2000l) observations, no obvious counterpart 
is visible that fully matches the morphology of this source. Here, 
we examine the possibility that HESS J1745-303 is (at least in 
part) associated with several counterpart candidates seen in other 
wavebands, most notably the interaction of a nearby supernova 
remnant with a molecular cloud. Additionally, we present an 
analysis of recent X-ray data from XMM-Newton covering the 
central part of this object. 



2. Technique 

2.1. The H.E.S.S. Instrument 

H.E.S.S. (the High Energy Stereoscopic System) is an array of 
four atmospheric Cherenkov telescopes (ACTs) located in the 
Khomas highlands of Namibia at an altitude of 1800 m above 
sea-level. Each telescope consists of a 107 m 2 optical reflector 
made up of segmented spherical mirrors that f ocus light into 
a cam era of 960 photo-multiplier tube pixels dBernlohr et al.l 
12003k Using the imag i ng atmospher i c Cherenkov techn ique 
(e.g. Hillasl [l 9851. 1 1 9961: IWeekesI [19961 iDaum et all 1 1 9971) . the 
telescopes image the Cherenkov light emitted by the particles in 
extensive air showers from multiple viewpoints, and the energy 
and direction of the primary gamma ray can be reconstructed 
with an average energy resolution of ~16% (above an energy 
threshold of approx imately 150 GeV), and a spatial resolution 
of ~0.1° per event ( lAharonian et al.ll2006al) . The large field of 
view (~5°), and good off-axis sensitivity of the H.E.S.S. array 
make it well suited for studying extended sources and for scan- 
based observations, where the source position is not known a 
priori. 



2.2. Data and Analysis Technique 

Following the standard H.E.S.S. procedure, the data pre- 
sented here are processed with separate anal ysis and calibra- 
tion sc hemes: the H.E.S.S. standard analysis (lAharonian et al.l 
l2006al) . in which showers are reconstructed and hadronic 
backgr ound is rejec ted via the Hillas moment-analysis tech- 
nique (|Hillas 1996b : and the Model2D analysis described by 
IdeNauroisI d2005l) . which employs a semi-analytic model of the 



shower to characterize each image. It should be noted that both 
the techniques and simulations used in them are independent, 
providing a robust check of the analysis. Since results of both 
analyses agree within errors, the results presented here are from 
the standard analysis only. The separation of gamma-ray candi- 
dates from cosmic-ray-like events was made using both stan- 
dard cuts (optimized for a lower energy threshold) and hard 
cuts (optimized for bet ter background rejection) described in 
( lAharonian et alj|2006al) . but with a larger angular angular in- 
tegration radius to account for the extension of the source. The 
former were used for spectra, and the latter for producing the sky 
images, though both were checked for consistency. 

The data for HESS J1745-303 are not primarily comprised 
of dedicated observations of the object, but rather from scan- 
based observations of the region (taken at regular grid-points 
along the galactic plane) and from dedicated observations of the 
VHE source at the Galactic center (HESS J1745-290), which 
lies approximately 1 .4° away. Due to the wide range of pointings 
and since the gamma-ray acceptance across the field of view of 
H.E.S.S. falls off radially, the total exposure of the region around 
HESS J1745-303 is highly non-uniform; in particular there is a 
strong gradient toward the Galactic center, which can lead to in- 
creased systematic errors in background subtraction. To mitigate 
this situation, several different background selection techniques 
were employed in this analysis. 

For the generation of two-dimensional images, both the field- 
of-view background method, where the background for each ob- 
servation is determined from a one-dimensional model of the 
radial acceptance (taken from observations with no significant 
emission in the field of view), and the ring background method, 
where the background at each point on the sky is calculated 
from an annulus surroundin g it (again, with sources excluded), 
are used ( jBerge et~aTll2007l) . Since the fie Id-of-view background 
method is more sensitive to gradients in exposure, the images 
presented here employ the ring background method, though both 
methods are checked for reasonable (within ~5%) consistency. 

For the spectral analysis background estim ation, we use 
the reflected region technique dBerge et al.ll2.007l) . where back- 
ground events are selected from circular off-source regions 
within the field of view. These regions are chosen with the 
same angular size and offset from the observation center position 
as the on-source integration region, ensuring their acceptance- 
corrected exposure is approximately equal to that of the on- 
source region. This technique is less suited to the generation of 
images, but for a known source position provides an estimate of 
the background that is independent of radial acceptance models. 

To prevent contamination of the background from the diffuse 
gamma-ray emission and to avoid including emission from the 
various other sources near the Galactic center, regions around 
Sagittarius A*, around G 0.9+0.1, as well as within ±0.8° of the 
Galactic plane, were not used for background estimation. 

The statistical significances in both the images and the spec- 
tral analysis are calculated from the measured number of on- 
and off-source (backgro und) events follo wing the likelihood ra- 
tio procedure outlined in lLi & Mai d!983l) . 

Spectra are generated fol lowing the methods described by 
both lAharo nian et al. (2006a), in the case of the standard anal- 
ysis, and iPiron et all (1200 ll) for the Model2D analysis. The on- 
source integration radius used for the generation of spectra is 
chosen (unless otherwise stated) to fully enclose the source 
(based on the radial profile of the signal and background), mak- 
ing no assumptions on the details of the source morphology and 
thus providing a less-biased result at the expense of a lower 
signal-to-noise ratio. Since the data set spans several years, dur- 
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H.E.S.S. 
(2006) 


H.E.S.S. 
(full) 


XMM-Newton 


Exposure Time (h) 


35 


79 


4 


Angular Res. (arcmin) 


6 


6 


0.1 


Field-of-view (°) 


5.0 


5.0 


0.5 


Energy Range 


0.35-1.5 TeV 


0.35-30 TeV 


0.5-10 keV 


Energy Res. 


16% 


16% 


20-50% 



HESS J1745-303 



Table 1. Chara cteristics of the observat ions of HESS J1745-303 
with H.E.S.S. dAharonian et al.ll2006al) and XMM-Newton. The 
"full" dataset is what is presented here, while the 2006 subset 
was presented bv lAharonian et all d2006dl) . The energy range for 
the H.E.S.S. case is the range used to fit the statistically signifi- 
cant spectral points; the full sensitivity of the detector extends to 
approximately 100 TeV. 



ing which the gain and optical efficiency of the telescopes has 
changed, muon images are us ed to calibrate the energ y estimate 
of each gamma-ray candidate dAharonian et al ]|2006al) . The sys- 
tematic error on the fluxes given here is estimated from simu- 
lated data to be +20% while the photon index has a typical sys- 
tematic error of +0.2. 



3. VHE Results 

With an exposure time of 35 hours, th e origina l H.E.S.S. 
data set presented by lAharonian et al.l d2006dl) revealed 
HESS J1745-303 at a pre-trials significance level of 6.3 stan- 
dard deviations (for an integration radius of 0.224°, which is 
the H.E.S.S. standard for blind source searches). With the in- 
creased exposure (now 79 hours) coming primarily from re- 
observations of the Galactic center source HESS J1745-290, 
HESS J1745-303 is now seen at well above the detection thresh- 
old, with a pre-trials significance level of 12 standard deviations 
{(f). 

Figure Q] shows an image of gamma-ray excess counts cover- 
ing HESS J1745-303 with significance contours overlaid. The 
image is smoothed with a Gaussian kernel with standard devia- 
tion 0.07° (chosen to reveal morphological features while main- 
taining good statistics), and the significance contours are gen- 
erated with an oversampling radius of 0.12°, matched to the 
RMS of the Gaussian smoothing to provide a visual impression 
of significant features. The source centroid is determined by fit- 
ting an elongated two-dimensional Gaussian convolved with the 
H.E.S.S. point-spread function to the un-smoothed images. Due 
to the non-Gaussian morphology, this only gives a rough cen- 
troid of the emission. 

For the spectral analysis, the integration region was taken 
to include the entire source (with a radius of 0.4° centered on 
17 h 45 m 2. 10 s , -30°22' 14.00", J2000 coordinates). For this larger 
integration region, we find a total significance of 10.2<x, with 
2030 total excess counts. The spectrum, shown in Figure [2] is 
well fit by a power law: dN/dE = Nq(E /TeV)~ r with photon 
index r = 2.71 + 0.11 sta t ± 0.2 sys , and a differential flux normal- 
ization N of (2.84 + 0.23 stat + 0.28 sys ) x lO-^TeV -1 cnrV 1 . 
This corresponds to an integral flux above the peak energy of 
VHE events (~1 TeV) of F(l - lOTeV) = (1.63 + 0.16) x 
10~ 12 cirT 2 s _1 . Figure|2]shows the resulting spectral energy dis- 
tribution. 

Figure [3] compares t he new result to the data presented in 
lAharonian etal.1 d2006dl) . Below 700 GeV, the new result no- 
ticeably differs from the former spectrum, which was made 
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Fig. 1. Top: A VHE gamma-ray (excess count) image of 
HESS J1745-303 with the positions of possible counterpart can- 
didates overlaid for reference. The color scale set such that the 
blue/red transition occurs at the ~ 3cr (pre-trials) significance 
level. The 4<x to 7<x statistical significance contours are shown 
in black. The thin white circles represent the integration regions 
A, B, and C discussed in the text. The dashed circle is the 95% 
error circle for the location of 3EG J1744— 301 1. The signif- 
icant excess seen to the north is th e tail end of the Galacti c 
Ridge diffuse emission discussed by lAharonian et alJ d2006d) . 
The Galactic plane is marked with a dotted line. Bottom: The 
H.E.S.S. s ignificance contour s (blue) overlaid on a VLA ra- 
dio image dLaRosaelall2000|l with overlaid ROSAT hard-band 
contours (green. fVoges et al1l2000l) . Stars show the positions of 
OH masers, and the XMM-Newton field of view (Figure HJ is 
drawn as a circle. 



with nearly a factor of 10 less photons (211 excess counts). 
Quantitatively, the new spectrum is 2.5cr softer than the old 
photon index of F = 1.82 + 0.29 sta t + 0.2 sys . Thus, the in- 
tegral flux above 200 GeV is smaller than the old value of 
(11.2 ± 4.0 stat + 3.4 sys ) x 10~ 12 cirrV 1 . 

A re-analysis of only the data used in the original publication 
with our updated procedures gives a result that is in agreement 
with the one presented here; and a re-analysis of all other sources 



4 



The H.E.S.S. Collaboration: HESS J1745-303 




1U 1U 

Energy (TeV) 



Fig. 2. The spectral energy distribution of VHE gamma rays 
from HESS J1745-303 (solid points) and from only region A 
(open circles). The data are fit by a power law of the form 
dN/dE = N Q (E/TeV)- T see Section E The 3EG J1744-3011 
flux is plotted as a bow-tie for comparison. 
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Fig. 3. Comparison of the spectra of HESS J 1745 -303 presented 
he re (filled circles with line fit) with the previous result presented 
by lAharonian et all (l2006dl) (open squares). 



presented by lAharonian et all d2006dh using the same techniques 
as in this paper gives results that are consistent with the original 
publication. Therefore, we have no indication that there is a gen- 
eral systematic error involved in either the old or the new analy- 
sis, but rather most likely a systematic error that exclusively (or 
predominantly) affects HESS J1745-303. 



Region 


F(l-lOTeV) 


% total 


r 




xl(T 12 cm- 2 s" 1 






Full 


1.63 ±0.16 


100 


2.71 ±0.11 


A 


0.25 ± 0.04 


15 + 3 


2.67 ±0.14 


B 


0.11 ±0.02 


6±2 


2.93 ±0.21 


C 


0.14 ±0.03 


8±2 


2.86 ± 0.27 



Table 2. Integral fluxes of the three test regions A, B, and C 
(shown in Figure [1), compared with the total integral flux from 
the full source. The photon index T is derived from a fit to the 
spectrum for each region. Note that region A is chosen to corre- 
spond with the molecular cloud position and radius described in 
section |5~T1 



When analyzing an extended source that has both low- 
surface-brightness and is located in a region of extremely uneven 
exposure, small uncertainties, e.g. in the acceptance-correction 
or background subtraction, may become significant. Such a 
problem only strongly affects low-surface-brightness sources an- 
alyzed with large integration radii. For most H.E.S.S. sources, 
wobble-mode observations are preformed that deliberately con- 
strain systematic background errors to a minimum; but in the 
case of HESS J1745-303, the data set is dominated by obser- 
vations of the Galactic center and the Galactic plane, therefore 
a uniform acceptance for all background control regions and 
the source region was probably not achieved in the old analy- 
sis to the required level. In contrast to the previous analysis, the 
large increase in observation time has provided a spectral-quality 
data set that is well above the detection threshold, over a much 
larger energy range, and that includes improvements in proce- 
dures used for the analysis of weak, extended sources. Here, 
we do not assume a particular source morphology and include 
time-dependent optical efficiency corrections to the energy, more 
exclusion regions for background subtraction (for sources that 
were subsequently discovered nearby), and software and lookup- 
table improvements that have reduced systematic errors due to 
uneven exposure. 

Since the emission appears to have a complicated morphol- 
ogy with more than one peak in the excess image, the possibility 
that HESS J1745-303 is more than one source was explored. 
First, three emission peaks were determined, located at the posi- 
tions A, B, and C shown in FigureQ] Between each pair of peaks, 
a one-dimensional slice in the uncorrelated excess image (with 
a width of 0.1°) was made to determine the significance of the 
"dip" between them. In each case, the emission is no more than 
two standard deviations from a constant value across the slice. 
Furthermore, if the emission peaks are from multiple sources, 
one might expect to see spectral variability across the object, 
though energy dependent morphology may also arise from trans- 
port and/or energy-loss processes within a single source. To test 
this possibility, a spectral analysis was made at each of the test 
points with an integration radius of 0.14° for A and 0.1° for B 
and C (see Table |2j. The spectral indices at each position are 
consistent with each other within statistical errors, and also with 
the spectrum determined for the entire source region; therefore 
within the statistics of the observations, there is no strong evi- 
dence to support the multiple-origin hypothesis. 



4. XMM-Newton Observations of the region 

HESS J1745-303 was observed with the XMM-Newton X-ray 
satellite on September 18, 2006 for 30 ks in satellite revolu- 
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Right Ascension 



Fig. 4. Point sources detected within the XMM-Newton ex- 
posure (see Table [3] for detailed information). The H.E.S.S. 
4, 5,6cr significance contours from Figure Q] are overlaid in 
white. Apparently associated stars are labeled. 
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Fig. 5. Exposure corrected X-ray image smoothed with a 0.01° 
Gaussian filter. All point sources above 3cr and the borders of 
each detector chip (shown by the contour) have been excluded 
before smoothing. No evidence for diffuse emission is seen. 



tion 1241 (ObsID 0406580201). All X-ray instruments (EPIC 
MOS 1, MOS 2, and PN) were operated in full-frame mode 
and a medium filter was applied to screen out bright optical 
and UV sources. The calibration, data reduction and analy- 
sis made use of the XMM-Newton Science Analysis Software 
(SAS), version 7.0, together with the Extended S ource Analysis 
Softw are package (XMM-ESAS), version 1.0 (ISnowden et al.l 
12004 . Following standard data reduction and calibration pro- 
cedures, the data set was cleaned from temporally occurring 



background caused by soft proton flares. The resulting obser- 
vation time amounts to 14.4 ks of useful data. Figure |4] shows 
an adaptively smoothed count map of the region surrounding 
HESS Jl 745-303 using events above 0.5 keV detected with 
either the MOS 1 or MOS 2 detector of XMM-Newton. The 
white contours indicate the H.E.S.S. VHE gamma-ray signifi- 
cance contours at 4<x, 5cr and 6<x. Six X-ray sources are appar- 
ent (labeled in red), as determined by the standard XMM-Newton 
source detection algorithm emldetect in the energy band from 
0.5-10 keV, as well as in sub-intervals from 0.5-2 keV, 2.0- 
4.5 keV and from 4.5-10 keV. Table [3] summarizes the sources 
detected above 0.5 keV. For all these sources the algorithm also 
attempts to determine a source extension by fitting a Gaussian 
model to the data. All six sources were found to be consistent 
with a point-source. Several of these X-ray sources coincide 
with stars known in the optical (shown in green in Figure [?}. 
Sources 1, 2, and 5 were found to be positionally coincident 
with HST optical guide stars as shown in Figure |4] source 3 
was found to coincide with an M3-star (RHI84 10-499) and 
for sources 4 and 6 a catalog search did not yield any obvious 
counterparts in other wavebands. As all of the detected X-ray 
sources are point-like and rather faint with fluxes around or be- 
low 10~ 14 erg crrT 2 s , assuming a similar energy flux in X-rays 
and gamma rays, it seems unlikely that any of these sources 
are connected to the bright extended VHE gamma-ray source 
HESS J1745-303, which is extended and has an energy flux of 
(2.00 + 0.18) x lO" 11 ergcirrV 1 above 200 GeV. In addition 
to the search for point-sources an analysis sensitive to diffuse 
X-ray emission has been performed. To that end, sources de- 
tected above a significance level of 3cr were excluded from the 
raw counts map and the exposure map. The maps were then 
smoothed with a Gaussian of width 0.01° and the ratio taken to 
produce the resulting smoothed, exposure-corrected counts map 
as shown in Figure [5] Diffuse emission at a level similar to or 
above the level of the detected point-sources should show up 
in this method. However, no sign of such a diffuse emission is 
detected in the whole field of view and we derive a 99% con- 
fidence limit on the flux level of the diffuse emission in region 
A of 4.5 x 10~ 13 erg cirT 2 s , using the rest of the field-of-view 
to determine the background, and 7.1 x 10~ 13 erg crrT 2 s using 
only a strip to the west of region A as background. 



0.05 5. Possible associations 



To look for possible associations, standard catalogs of sources 
thought to be associated with VHE em ission were searched, in - 
cluding high-spin-do wn flujfl pulsars ([Manchester et al. 2005 ), 

"1, 



SNRs dGreenl 12004 . Wo lf-Rayet stars (Ivan der Hucht 



2001 



high-mass X-ray binaries dLiu et al.ll2006l) . INTEGRAL sources 
dBird et al.l 120071) . and HII regions. Ad ditionally , public sur- 
vey d ata from ROSAT dVoges et alj|2000l), A SCA dTanaka et al.l 
119941) . and the VLA (LaRosa et al.l 120001) were searched for 
possible un-cataloged counterparts. The most likely candidates 
found in this search are discussed here. 



5.1. G359. 1-0.5 and Molecular Clouds 



G359. 1-0.5 was identified as an SNR bv lDownes eTail d 19791) 
using the Westerbork Synthesis Radio Telescope and the Bonn 



1 High spin-down flux (E/D 2 ) pulsars are thought to be likely VHE 
emission candidates if the conversio n efficiency from spi n-down power 
to VHE emission is around 1% (e.g. lCarrigan et ai1l2007h . 
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Fig. 6. Velocity-in tegrated (-100 to -60kms _1 ) map of 12 CO 
dBitran et alJI 19971) emission showing the molecular cloud coin- 
cident with the shell of G359. 1-0.5, with overlaid H.E.S.S. 4<x 
to 7cr significance contours. 



100-m telescope. VLA observations dUchida et alj|l992al) con- 
firmed the presence of the non-thermal shell and HI absorption 
showed that its probable location is within a few hundred par- 
sees of the Galactic center. At an estimated distance of 7.6 kpc, 
the 12' radius of the remnant corresponds to 26.5 pc, with an 
estimated age of > 10 4 yr (a middle aged SNR). Observations 
of the 12 C O emission line (J=l — >0) with the Bell Laboratories 
telescope dUchida et al.l 1 1 992bh revealed a dense ring of mat- 
ter surrounding the shell. The radial velocity dispersion of this 
super-shell between -60 and -190kms 1 agrees well with a lo- 
cation of the remnant close to the Galactic center. This super- 
shell could have been produced by the combined stellar winds 
of ~200 O-type stars concentrated in this region , which is ver y 
probably the birth place of the remnant dUchida et al.l 1 1992bb . 
X-ray observations of the remnant with ASCA showed no evi- 
dence for a shell in the energy range 3.2-10.0keV, detecting only 
therm al diffuse emission from the central region dBamba et alj 

120011 

lUchida et~ai1dl992bl) observed a correlation in intensity be- 
tween the non-thermal radio emission of the remnant and the CO 
emission and suggested an association of the two. There are no 
sharp gradients observed in the velocity distribution, the absence 
of which could be explained by the presence of magnetic precur- 
sors or by a previous acceleration of the shell by stellar winds or 
old supernovae. Another indication of the interaction of the rem- 
nant with the surrounding medium is the presence of maser emis- 
sion s pots near the edge of the shell, revealed by VLA observa- 
tions dYusef-Zadeh et al.ll 19951) . OH masers at 1720 MHz are be- 
lieved to be produced by collisional pumping beh ind shocks an d 
are therefore a good indicator of shocked clouds dElitzurll 19761) . 
Although the mean velocity of the masers at around -5 kms~' 
is significantly shifted from the velocity range of the 12 CO shell 
([-60,-190] kms 1 ), a random c oincidence of the shell wit h the 
maser positions seems unlikely. lYusef-Zadeh et aL I (119951) dis- 
cussed this velocity discrepancy: the two velocities could be rec- 



onciled if the shocked cloud was carried away by the shock it- 
self, which could have redirected the radial velocity component. 
The specific conditions required to allow a population inversion 
(and thus maser emission) may also explain this discrepancy. 
Moreover, the association between the OH masers, the SNR and 
the 12 CO cloud is supported by the fact that the maser distribu- 
tion shows a good correlation with the 12 CO emission maximum 
and the non-thermal radio emission from the western part of the 
shell. 

We have used 12 CO (J=l->0) data from the Cerro Tololo 
Inter-American Observatory (Chile) to ma p the matter distri - 
bution in the western part of the remnant dBitran et al.l 1 1 997b . 
Figure [6] is the 12 CO map integrated in velocity between 
-100 kms 1 and -60kms _1 . This velocity range corresponds 
to the western part of the super-shell surrounding the remnant 
and is potentially associated with the remnant. Two overlap- 
ping clouds are present in this region. The cloud shown in 
Figure [6] contains a component that is partially coincident with 
HESS J1745-303. 

Due to the relatively high magnetic field, (0.2-0.6 mG, 
as measured from the Zeeman splitting of maser lines 
dRobinson et al.fl996l) . the cooling time of ultra-relativistic elec- 
trons would be much shorter than the age of the remnant, and 
therefore an electron production scenario seems unlikely as an 
explanation of the VHE emission. A possible explanation is 
that a part of the VHE gamma-ray source comes instead from 
hadronic cosmic-ray interactions in this cloud, producing neutral 
pions that decay into two gamma rays. Assuming a value of the 
ratio X = N Hl /W C o of l.SxlO^cm^K 1 km" 1 s 1 dDame et al.l 
1200 ll) . we estimate the H2 mass in this cloud to be 5 x 10 4 M o 
with a density of 5 x 10 3 cirT 3 , assuming spherical symmetry. 
The value of the ratio X used is a mean galactic value which 
may not be relevant for the inner part of the galaxy, therefore the 
cloud mass may be overestimated. A spectral analysis of the part 
of HESS J1745-303 coincident with the cloud (region A) indi- 
cates that approximately 15% of the observed flux comes from 
this region (see Table|2]i. 

Since hadronic interactions lead to the production of gamma 
rays with energies typically a factor of 10 below the primary en- 
ergy, the energy of interacting protons (W) required to generate 
the observed flux of VHE y-rays between 300 GeV and 40 TeV 
can be estimated in the corresponding energy range of approx- 
imately 3-400 TeV to be: W = t pv ^ n o ■ L y , where ? pp ^ is the 
characteristic cooling time of protons through the 7r° production 
channel and L y is the gamma-ray lumin osity between 300 GeV 
and 40 TeV dAharonian & Atovanl 19961 ). Assuming that the pro- 
ton energy distribution follows a power law with the same index 
as the gamma rays over the full relativistic range, we can ex- 
trapolate this distribution down to 1 GeV. Assuming also that the 
accelerated proton density is uniform in the whole remnant, it 
corresponds to a fraction of ~32% of the mechanical explosion 
energy of 10 51 erg of the remnant. This estimate suffers from 
large uncertainties, mainly from the cloud mass estimation, the 
fraction of HESS J1745-303 involved in this association, and the 
explosion energy. However, it is interesting that we obtain an es- 
timate that is comparable to the theoretical espectation of -10%. 

5.2. 3EGJ1 744-3011 

Most of the emission associated with HESS J1745-303 lies 
within the 95% error-circle of the unidentified EGRET source 
3EG J1744— 301 1 (marked with a dotted circle in Figure [T), 
which has an integral flux in the 100 MeV-10 GeV energy 
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band of (63.9 + 7.1) x 10~ 8 cm -2 s~' with a photon index of 
Tegret = 2.17 + 0.08 ( lHartman et al.lll999l) . Extrapolating this 
flux to the VHE range, we find an expected integral flux in the 
range 1-10 TeV of x 10~" cirT 2 s , which is higher than 
the integral flux observed by H.E.S.S. over the same range (see 
Figure |2j. Fitting a power law with an exponential cutoff to the 
H.E.S.S. spectrum combined with the EGRET flux point, we find 
a cutoff energy of ~(3.0 + 1.5) TeV. 

A study of the long -term variability of EGRET sources made 
bv lTorres etail d200ll) shows that 3 EG J1744-3011 is variable 
(on year timescales), with an average statistical index of vari- 
ability over 3cr higher than that expected from pulsars (which 
are considered a non-variable source used as a reference for sys- 
tematic variability of the instrument). If the EGRET source is 
truly variable, it is unlikely to be associated with the extended 
emission seen in HESS J1745-303, where no variability would 
be expected. 



5.3. PSR B1 742-30, PSR J 1747-2958, SLX 1 744-299 

Two cataloged pulsars lie within or near HESS J1745-303: 
PSR B1742-30 and PSR J1747-2958 (see FigureQ). Energetic 
pulsars driving pulsar wind nebulae are known to produce VHE 
emission t hat may be asymmetric o r offset from the pulsar po- 
sition (e.g. lAharonian et aLll2006b 1l5 l2005d llbh. In several cases, 
the PWN candidate has been first identified in the VHE energy 
rang e and subsequently confirmed with X-ra y measurements 
(e.g. lAharonian et al.ll2007t iHinton et alj|2007l) . 

PSR B 1742-30 (also known as PSR J 1745-3040) is a rather 
old pulsar (546 kyr) located near the southern edge of the re- 
gion of significant emission in HESS J1745 -303, and has a 
spin-d own flux E/D 2 — 2 x 10 33 erg s _1 kpc 2 dManchester et alj 
120051) . requiring a conversion efficiency from rotational kinetic 
energy to gamma-ray emission of approximately 32% to pro- 
duce the entire VHE emission. Such high apparent efficiencies 
are possible if the spin-down flux was much higher in the past 
and the particle cooling times are comparable to or much shorter 
than the pulsar age. If PSR B1742-30 powers only the fraction 
of HESS J1745-303 enclosed by region B, then the required 
conversion efficiency would be only 2%, which is not unreason- 
able compared with other known VHE PWNe. 

PSR J1747-2958 is located approximately a half-degree east 
of HESS J1745-303 and is associated with the bright X-ray 
and radio feature G3 59.23-0.82, also known as "the Mouse" 
dGaensler et al.l 120041) (seen in the radio image in Figure Q]), 
which is believed to be a bow-shock PWN with a trailing tail 
caused by the reverse termination shock . Given the proper ve - 
locity, estimated to be ~ 600 km s -1 dGaensleret al. 2004), 
distan ce D 2.08 kpc and age (25.5 kyr) dManchester et alj 
120051) . one can extrapolate that the pulsar would have moved 
approximately 0.43° from its original position, placing it close 
to HESS J1745-303 if the direction of motion is along the 
tail of the "Mouse" (see dotted line in Figure Q] Bottom). The 
relatively high spin-down flux of PSR J1747-2958 (E/D 2 = 
4 x 10 35 erg s _1 kpc -2 ) would imply a 0.2% conversion efficiency 
to explain the entire VHE emission, or 0.02% for only region 
A. In this case, the PWN would have to be extremely offset and 
asymmetric. 

Located near PSR J1 747-2958 are the ultr a-compact X-ray 
binary SLX 1744- 299 din 't Zand et al.l 120071) . and the X-ray 
burster 1744 - 300 dSkinner et aljl990b . which are not generally 
expected to produce offset or extended emission and are thus not 
considered probable counterparts to HESS J1745-303. 



6. Discussion 

Due to the positional coincidence and plausible energetics, at 
least part of the emission of HESS J1745-303 (in the region 
labeled A in Figure []]) may well be associated with the inter- 
action of the shell of G359. 1-0.5 with a molecular cloud. This 
scenario fits particularly well within the context of theoretical 
predictions for VHE gamma-ray emission from SNRs embed- 
ded in dense media. iGabici & Aharonianl ([2007) show that for 
a SNR of approximately the same age as G359. 1-0.5 in this 
context, one would expect significant TeV emission that peaks 
around 1 TeV. Assuming a similar supernova energy output of 
10 51 erg and scaling their theoretical model for a cloud at 30 pc 
away from an SNR shell by M c \/D 2 , where M c \ is the mass 
of the cloud near G359. 1-0.5 and D is its distance, we find 
that the predicted flux at 1 TeV of 2.8 x 10~ 12 TeV" 1 crrT 2 s" 1 
matches well with the HESS J1745-303 emission around region 
A, F A (1 TeV) = (4.4 + 0.5) x 10~ 12 TeV" 1 cirT 2 s" 1 , as does the 
soft T ^ 3.0 spectral index. 

However, since the VHE emission extends beyond region A 
(and the dense target material does not), the SNR/molecular- 
cloud scenario is not sufficient for describing the entire VHE 
source. Since there is as yet no statistically significant sep- 
aration between the various parts of the emission region, an 
explanation for the entire source remains complicated and 
the possibility of source confusion still remains. In particular, 
PSR B 1742-30 is energetic enough to power a PWN in the part 
of HESS J1745-303 surrounding it (region B), and though an 
unlikely candidate due to its significant offset, PSR J1747-2958 
is in principle powerful enough to power the entire VHE source. 
The lack of significant spectral variability across the emission 
region further complicates the identification of counterparts. It 
is important to note that many VHE sources and source classes 
have spectral indices in the range 2.1-2.7, therefore given the 
statistics, it would be likely impossible to disentangle a superpo- 
sition of three randomly chosen VHE sources from their spectra 
alone. 

Furthermore, the lack of a significant extended source in the 
XMM-Newton data is not without prec edent — several o t her TeV 
sources, such as HESS-J130 3-631 dAharonian et alj 120053). 
other unidentified TeV sources dAharonian et al.ll2008l) . have so 
far no identified X-ray counterparts The lack of non-thermal X- 
ray or radio emission combined with the relatively high magnetic 
field around G359. 1-0.5 further supports the hadronic scenario 
for VHE gamma-ray production for the region associated with 
the interaction of the SNR shock with the target material in the 
molecular cloud. However, even in a hadronic scenario, some 
longer-wavelength emission would be expected due to secondary 
electrons produced in the interactions. 

An association of 3EG J1744-3011 with part or all of 
HES S J1745-303 is al so plausible from an energetic standpoint 
(e.g. iFunk et alj 120071) . however since the size of the EGRET 
error circle is larger than the VHE emission, the position may 
not correspond with the VHE source and further lo calization 
is not currently possible. lAharonian & Atovanl £l996) show that 
the observed gamma-ray flux from a hadronic source is propor- 
tional to E y P+S , where Y p is the proton index at the source and 
6 is the index of the diffusion coefficient (typically 0.3 - 0.6), 
allowing spectra that are quite soft in the TeV energy range, 
and w hich may be diffe rent from the slopes in other energy 
bands dTorres et alj|2003l) . Therefore for an E p 2 source spec- 
trum, it is possible to reproduce both the hard (F = 2.17) EGRET 
spectrum and the softer F = 2.71 + 0.11 sta t ± 0.2 sys TeV spec- 
trum. If we assume 3EG J 1744-301 1 is associated only with the 
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SNR/molecular-cloud interaction, then the observed flux is sig- 
nifican tly higher than that predicted by the iGabici & Aharonianl 
(120071) model, however given the uncertainties in position, the 
EGRET flux may also contain contributions from the SNR shell 
itself (where gamma rays may be produced by e.g. inverse- 
Compton scattering of lower energy photons by high-energy 
electrons accelerated in the shock), or from other sources in the 
region, which could explain this. 

Further multi-wavelength observations in X-rays and GeV 
gamma rays (e.g. from observatories such as XMM, Suzaku and 
GLAST) of the entire region spanned by HESS J1745-303 as 
well as deeper VHE exposures will be needed to disentangle the 
emission possibilities and to find more definitive counterparts. 
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Table 3. Sources detected using the detection algorithm emldetect. The parameters given here are for the energy range between 
0.5 and 10 keV. Col. 2: Name recommended by the XMM-Newton SOC and the IAU for source detections. Col. 3 and 4: J2000.0 
coordinates. Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arc-minutes, and arc- 
seconds. Col. 5: Error on the source position in arc-seconds. Col. 6: Number of counts in EMOS1 and EMOS2 within a 10" 
integration region using events above 0.5 keV. Col. 7: Statistical significance of the detection derived with emldetect. 8: Source flux 
above 0.5 keV in 10~ 14 erg cirT 2 s _1 . 



Id 


XMMUJ 


RA 2 ooo 


D ec 2000 


2-D Error 


Counts 


Significance 


Flux-10- 14 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1 


174434.5-301522 


17 h 44 m 34.47 s 


-30°15'21.9" 


0.6 


280 


14.3 


10.4 ± 0.7 


2 


174441.3-301648 


17 h 44 m 41.27 s 


-30° 16' 47.8" 


0.6 


209 


12.0 


8.4 ±0.7 


3 


174437.6-301812 


17 h 44 m 37.55 s 


-30° 18' 12.3" 


1.3 


51 


5.2 


2.1 ±0.4 


4 


174425.3-302236 


17 h 44 m 25.30 s 


-30°22'36.1" 


1.3 


70 


7.0 


3.2 ±0.4 


5 


174458.6-301743 


17 h 44 m 58.59 s 


-30°17'43.2" 


0.9 


95 


7.5 


5.2 ±0.7 


6 


174351.0-301709 


17 h 43 m 5 1.03 s 


-30°17'09.0" 


1.0 


76 


6.9 


4.4 ±0.6 



